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ABSTRACT 



We consider a Hartmann layer, stationary flow of a viscose and resistive fluid be- 
tween two plates with superimposed transverse magnetic field, in the limit of gyrotropic 
plasma, when viscosity across the field is strongly suppressed. For zero cross-field vis- 
cosity, the problem is not well posed, since viscosity then vanishes on the boundaries and 
in the middle of the layer, where there is no longitudinal field. An additional arbitrarily 
small isotropic viscosity allows one to find magnetic field and velocity profiles which 
are independent of this viscosity floor and different from flows with isotropic viscosity. 
Velocity sharply rises in a thin boundary layer, which thickness depends both on the 
Hartmann number and on the Lundquist number of the flow. 

The implication of the work is that, in simulating ICM dynamics it is imperative to 
use numerical schemes which take into account anisotropic viscosity. Although magnetic 
fields are dynamically subdominant in the ICM they do determine its the dissipative 
properties, stability of embedded structures and the transition to turbulence. 

Subject headings: galaxies: clusters: general 



1. Introduction 



Gasdynamical interactions of magnetized flows in the cores of clusters of galaxies play an 
important role in formation of the observed morphological structures, e.g. through expansion of 
AGN blown bubble into ICM medium and resulting plasma heating, and interaction of two gas 
components in merging clusters (e.g. Markevitch Vikhlinin||2007 ). Majority of theoretical work 
on these topics has been numerical, mostly using the existing fluid and MHD codes, like ZEUS. Un- 
fortunately, simple hydrodynamic models "face multiple failures" ( Reynolds et al.||2005 ) for cluster 
cores. Perhaps, the most evident example is the Raleigh- Taylor and Kelvin-Helmholtz instabilities 
of AGN blown bubbles, which disrupt the bubbles approximately on one rise time (e.g 



Kaiser 



et al. 2005|. On the other hand, artificial fiddling with viscosity (which is usually parametrized 



with respect to Braginsky (1965) value; this procedure is not justified in ICM (Schekochihin et al. 



2005 Lyutikov 2007 1), shows that "modest" changes of shear viscosity lead to qualitatively differ- 



ent results (e.g. Reynolds et al. 2005; Sijacki & Springel 2006, increased viscosity make the ICM 



plasma gel-like and quenches the instability). 
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One of principal reasons, perhaps, for the failure of simple MHD codes is that they use isotropic 
viscosity, while ICM plasma is strongly gyrotropic, in a sense that it is weakly collisional, tl <C A < 
L (tl is ion Larmor radius, A is mean free path and L is a typical size of the system), and weakly 
magnetized, kinetic pressure p is much larger than magnetic pressure, (3 = 8irp/B 2 S> 1 (B is a 
typical magnetic field ). In ICM the ion Larmor radius ri ~ 10 8-9 cm is many orders of magnitude 
smaller than the mean free path A ~ 10 22 ~ 23 cm and the size L ~ 10 24 cm, while /3 ~ 100 (Carilli 



k Taylor 2002). 



In a strongly gyrotropic plasma the local transport properties, primarily viscosity and conduc- 



tivity, become anisotropic (Braginsky 1965). The effects of anisotropic viscosity and conductivity 



are expected to change substantially the results of ICM simulations. As a simple example, note 



that magnetic field drapes around the contact surface separating the two interacting media (Lyu 



tikov 2007). In a strongly gyrotropic plasma the shear viscosity inside the draping layer, with a 



flow along magnetic field lines, then becomes 0. This runs contrary to the idea that high viscosity 



may provide stabilization Reynolds et al. ( 2005 ) . (Note that draping itself can provide stabilization 



against KH instability, in a way similar to the effect of a thin oil layer on water waves ( Dursi|2007 )). 



Understanding the basic properties of strongly gyrotropic plasma is imperative for further 
progress, especially for parametrization of the 'sub-grid' physics in large numerical simulations. In 



this Letter we adapt one of the basic solutions of MHD, a Hartmann flow (e.g. Landau & Lifshitz 



1984), to the case of anisotropic viscosity. 



2. Hartmann flow with anisotropic viscosity 



Consider a one-dimensional (along x direction) flow of a weakly collisional plasma between two 
plates located at z = dba with a superimposed external magnetic field (generally, oblique) . This 
is meant to represent a boundary layer during interaction of two plasma flows in the ICM. In the 



Chew-Goldberger-Low (zero Larmor radius and neglecting heat fluxes) approximation (Chew et al. 



1956) the equations of resistive plasma flow read (Kulsrud 2005): 



^ = -V (P x + B 2 /2) + V (SS(P|| - Pi - B 2 

PH - p ± = 3 m dt lnB = 3i] (bbW 
<9 t B = curlv x B + rj r ■ AB 



(1) 



where t]q is the first Braginsky coefficient (Braginsky 1965), resistivity r\ r is a tensor and b = B/f? 
is a unit vector along magnetic field . We also absorbed a factor y/4ir into definition of magnetic 
field, 



Assuming that all quantities are independent of x and y, from divB = we find B z = const = 
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Bq. Introducing B x = BqJ(z), v x = v, the x component of Euler equation gives 

Blf + 3nod g ((YT^p) " ^ = ^ 

d x P± = AP/L is constant pressure gradient driving the flow, AP is a drop in pressure over the 
length L. Eq @ can be integrated once: 

Blf - A ^(z-z ) + 3vo = (3) 

Integration constant zq is a point where / = 0. 

Assuming that the y component of magnetic field vanishes, so that current flows across mag- 
netic field j = j y e y , the resistivity equation gives: 



v 



' = ~V±f" (4) 



Here i]± is resistivity across magnetic field . Finally, eliminating v' from Eq. ^ using Eq. ^ , we 
get equation for /: 

A.P f" f 2 

B of ~ -jr-(z ~ z o) - 3wu ^ - py = (5) 

This is the main equation which determines the structure of the flow. In case of anisotropic viscosity 



it's a nonlinear equation. It differs from the case of isotropic viscosity {e.g. Landau Sz Lifshitz 



1984 §67) by a different viscosity term: 



'(1 + / 2 



r] eff = Sr]o n | (6) 



Renormalizing, z — > za, zq — » zqcl and introducing a = {B 2 / AP){L / a) , we find 

z °~ z + a { f ~^(iTJ 2 T) =0 (7) 

where ^ 

G V^VOVr 8 ^ 

we identify as Hartmann number. 

Boundary conditions are somewhat tricky in our case. From the continuity of the tangential 
magnetic field (assuming there is no surface current) and from the symmetry of the flow, it is 
required that / = /o at z = ±a and / = at the point zq, where /o is an imposed longitudinal 
magnetic field . Below we consider a case when there is no superimposed longitudinal magnetic 
field, /o = and z$ = 0. 

On the other hand, at surfaces where B x oc / = (at the boundaries and in the middle) there 
is no viscosity, see Eq. Q, so the usual conditions of v = at z = ±a need not to be satisfied. 
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Thus, in principle, the flow may slip along the boundaries and have a discontinuity in the middle. 
This will make the problem unsolvable as the order of the ODE would be higher than the number 
of boundary conditions. 

In fact, it is necessary to assume some floor for viscosity to get a physically meaningful solution 
ofEq. @ even if we just impose condition f = at z = =ta. Near points / = 0, Eq. ([7| reduces to 
/ = {z — zq)/ol. So that derivative of / at these points has a definite sign, given by the parameter 
a. Thus, / cannot be zero more than once. Since / = in middle of the layer, this would clearly 
violate the conditions that parallel magnetic field vanishes on the walls. 



Introducing resistivity floor f r in Eq. §7j by substitution f 2 — > f 2 + f r , we can integrate Eq. 
([7]) numerically, see Fig.Q. For sufficiently small values of this viscosity floor the final result is 
independent of its exact value, Fig. ([2]). 

The somewhat unphysical value of can be expressed in terms of bulk velocity in the middle 
of the layer vq: 

AP vqBI _ vorjo 
L ~ G 2 rj ± ~ a 2 

This gives a = = G 2 /L u , where we identified the ratio VQa/r]± with magnetic Reynolds, or 

Lundquist number L u . Parameter a plays an important role in determining dynamics of the layer. 
If viscosity is dominated by ion-ion collisions, we may estimate a ~ l/({3M s K n ), where (3 is plasma 
beta parameter (ratio of kinetic to magnetic pressure), M = Vq/c s is the flow Mach number (ratio 
of velocity to sound speed c s ), K n = A/a is Knudsen number, the ratio of the mean free path A to 
the characteristic length scale a. 



0) 



These solutions are quite different from the case of isotropic viscosity ( Landau fe Lifshitz| 1984 1 



and have a number of particular feature. On the one hand, for large G > 1, in the bulk of the 
flow the profile of / is linearly increasing, / ~ z/a (flat velocity profile), similar to the isotropic 
viscosity case. The parallel magnetic field drops back to within a narrow boundary layer. Let us 
estimate the thickness of this boundary layer. Near the boundary / — > 0, so that an approximate 
solution of Eq. ([7| is 

^_Z a /z\(l+Vl+4G2a2 )/2 



a a \CL' 

If we define the thickness of the boundary layer 5 when /' = 0, we find 

S ln(aG) , .L u 

^ ln(aG )-g 11 

a aG G^ 

where we assumed Ga S> 1. This expression can be compared with the case of isotropic viscosity, 
when 5i/a ~ 1/G. The ratio 5/5i ~ 1/a = L u /G 2 . Thus, parameter a = G 2 /L u measures the 
relative concentration of a magnetic field profile toward the wall; for a > 1 the boundary layer in 
case of anisotropic viscosity is narrower than in the isotropic case. 

In most applications the ratio L u /G 3 is very small. For example, for the typical parameters of 
intercluster medium, velocity ~ 1000 km/s and layer thickness ~ lkpc, plasma beta j3 = 100, we 
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Fig. 1. — Function / (longitudinal magnetic field) for a = 1 and different G = 1, 3.1, 4.5 (bottom 
to top) (larger values of G are difficult to resolve numerically). 
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z/a 



Fig. 2. — Effect of finite isotropic viscosity on the structure of magnetic field for f r = 
0.1, 0.005, 0.0001, 1CT 10 (bottom to top), a = G = 1. For small enough f r solutions are nearly 
independent of the exact value of f r . 

estimate L u = 10 27 , G = 10 12 , so that L u /G 3 ~ 10~ 9 , thus S < a. 

3. Conclusion 

In this letter we considered a basic plasma physics problem: Hartmann flow with anisotropic 
viscosity. We first argued that in the case when transverse viscosity is suppressed completely, the 
problem cannot be formulated in a physically meaningful way: there should be some small isotropic 
contribution to viscosity. We derived magnetic field and velocity profiles, which in the limit of small 
isotropic viscosity floor are independent of its exact value. These profiles are considerably different 
from the case of isotropic viscosity. The velocity gradients are much more concentrated close to 
the walls of the channel than in the case of isotropic viscosity. 

How important is the structure of a boundary layer for the overall structure of the flow? On 
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Fig. 3. — Comparison of analytical magnetic field profiles (multiplied by a for better graphical 
representation) for isotropic (Landau & Lifshitz 1984 / oc {{z/a) sinhG — sinh(zG/a)) /(coshG — 
1)) and approximate anisotropic (Eq. (|10[)) viscosity for G = 10 and different parameters a = 



0.1, 1 10. For a > 1, the boundary layer in case of anisotropic viscosity is narrower than in the 
isotropic case. 
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the one hand, in a laminar regime at low Reynolds numbers, the structure of the boundary layer is 
probably not important: in the boundary layer the relative velocity just drops to zero, according to 
some law, without affecting the overall structure of the flow. On the other hand, properties of the 
boundary layer determine its stability and, thus, transition to turbulence ( Landau Lifshitz|19 59). 



In ICM plasma the Reynolds numbers are in the range Re ~ 10 — 1000, while, typically, transition 
to turbulence occurs at Re ~ 100 — 1000. (Transition to turbulence occurs in the boundary layer). 
Thus, we expect that effects of anisotropic viscosity are likely to be very important for ICM plasma, 
especially for determining its transition to turbulence. 

Thus, in simulating ICM dynamics it is imperative to use appropriate numerical schemes which 
take into account anisotropic viscosity (and to a lesser extent, conductivity), like the ones that 
have been applied to accretion disks ( Sharma et al.|[2~007 ). Though magnetic fields are dynamically 



subdominant in ICM (plasma beta parameter is large), they do, in fact, determine its the dissipative 
properties, stability of embedded structures and transition to turbulence. 
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